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Abstract 

It  is  shown  that  for  four-transmitter  systems,  a  family  of  four-by-four  unit-rate  complex  quasi-orthogonal  space-time 
block  codes,  where  each  entry  equals  a  symbol  variable  up  to  a  change  of  sign  and/or  complex  conjugation,  can  be 
generated  from  any  two  independent  codes  via  elementary  operations.  The  two  independent  groups  of  codes  in  the  family 
generally  have  different  properties  of  diversity,  but  the  codes  in  each  group  have  the  same  diversity  provided  that  the 
differential  symbol  constellation  is  symmetric.  It  is  also  shown  that  for  four- transmitter  systems,  an  eight-by-four  unit-rate 
complex  linear  dispersion  space-time  block  code  can  be  constructed  by  using  Hurwitz-Radon  families  of  matrices  of  size 
eight  such  that  diversity  three  is  guaranteed  even  when  all  symbols  are  independently  selected  from  any  given  constellation. 
This  code  is  so  far  the  only  known  unit-rate  linear  dispersion  code  that  has  diversity  no  less  than  three  for  four  transmitters 
under  any  given  constellation. 

©  2008  Elsevier  B.V.  All  rights  reserved. 

Keywords:  Space-time  block  codes  (STBC);  Orthogonal  STBC;  Quasi-orthogonal  STBC;  Non-orthogonal  STBC;  Hurwitz  Radon 
families  of  matrices;  Diversity  analysis 


1.  Introduction 

Design  and  analysis  of  space-time  block  codes  (STBC)  for  multiple  transmitting  antennas  have  been  an 
active  field  of  research  since  the  work  by  Alamouti  [1]  and  that  by  Tarokh  et  al.  [2].  STBC  is  aimed  to  exploit 
the  channel  diversity  between  multiple  transmitters  and  multiple  receivers  to  improve  the  rate  of  reliable  data 
transmission  and/or  the  performance  of  bit  error  rate.  STBC  is  also  useful  for  cooperative  relays  in  wireless 
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mobile  networks  [3-6],  where  STBC  can  be  used  effectively  as  if  between  multiple  transmitters  and  a  single 
receiver. 

STBC  is  a  mapping  (applied  at  the  transmitters)  between  a  sequence  of  input  symbols  and  multiple 
sequences  of  output  symbols.  The  number  N  of  the  output  sequences  typically  corresponds  to  the  number  of 
transmitters.  The  ratio  of  the  length  T  of  the  output  sequences  over  the  length  S  of  the  input  sequence  is  called 
the  rate  of  the  STBC  (assuming  that  both  the  input  symbol  constellation  and  the  output  symbol  constellation 
have  the  same  dimension).  The  output  of  the  STBC  mapping  can  be  denoted  by  a  T  x  N  matrix  C(v)  where 
the  S  x  1  vector  v  represents  the  input  symbols.  Assume  that  the  channel  following  the  N  transmitters  is 
frequency  flat,  and  there  are  M  receiving  antennas  at  the  end  of  the  channel.  Then,  the  received  baseband 
signals  at  the  destination  over  a  time  interval  of  T  symbols  can  be  represented  by  the  T  x  M  matrix  Y: 

Y  =  C(x)H  +  W,  (1) 

where  H  is  an  N  x  M  channel  matrix  whose  entries  may  be  assumed  to  be  i.i.d.  complex  Gaussian  random 
variables  (Rayleigh  fading),  and  W  is  a  T  x  M  noise  matrix  whose  entries  may  also  be  assumed  to  be  i.i.d. 
complex  Gaussian  random  variables.  With  a  coherent  maximum  likelihood  decoder,  the  pairwise  error  rate 
(PER)  P(x  ->  x)  averaged  over  the  channel  fading  distributions  is  upper  bounded  as  follows  [2]: 

E „[P(x  v/j  (Es/4N0)-rM,  (2) 

where  Es  and  No/2  are,  respectively,  the  symbol  energy  and  the  variance  of  noise  per  dimension;  the  signal- 
to-noise  ratio  (SNR)  may  be  defined  as  the  ratio  of  Es  over  No/2;  r  is  the  minimal  rank  of  C(x  —  x)  over  all 
possible  distinct  pairs  of  the  symbol  sequences;  vj  (j=\,...,r)  are  the  non-zero  eigenvalues  of 
C(x  —  x)HC(x—  x).  To  reduce  the  PER  of  a  code,  one  must  increase  r  and  the  minimum  of  Ylj=\vj-  The 
value  of  r  is  called  the  diversity  of  the  code,  and  the  minimum  of  Ylj=\vj  determines  a  coding  gain.  Diversity 
and  coding  gain  are  among  the  key  measures  of  a  code. 

A  detailed  review  of  STBC  is  available  in  [7,8].  For  convenience,  we  will  also  refer  to  STBC  simply  as  codes. 
The  most  attractive  codes  are  perhaps  the  orthogonal  codes  [1],  which  allow  the  maximum  likelihood 
(optimal)  detection  to  be  performed  independently  on  each  of  the  individual  symbols.  But  the  unit-rate 
orthogonal  complex  codes  exist  only  for  two  transmitters  [9].  For  more  than  two  transmitters,  there  are  only 
fractional-rate  orthogonal  complex  codes  [9].  Upper  bounds  on  the  rate  of  orthogonal  complex  codes  are 
explored  in  [10].  There  are  also  quasi-orthogonal  codes  that  allow  the  maximum  likelihood  detection  to  be 
performed  independently  on  pairs  of  symbols  [1 1]  or  even  independently  on  each  symbol  as  shown  in  [12].  But 
the  quasi-orthogonal  code  given  in  [11]  does  not  have  a  full  diversity.  Various  improvements  of  quasi- 
orthogonal  codes  are  further  developed  in  [13-16].  In  [16],  it  is  shown  that  unit-rate  quasi-orthogonal  codes 
with  maximal  diversity  products  can  be  constructed  by  using  a  finite  information  symbol  set  on  square  and 
triangular  lattices.  There  are  also  codes  that  are  designed  to  maximize  an  orthogonality  measure  [17]. 
Numerous  other  codes  can  be  found  via  [18-24]  and  the  references  therein. 

The  purpose  of  this  paper  is  not  to  present  a  new  code  competing  against  existing  ones.  But  rather,  we  reveal 
a  structural  insight  into  a  class  of  linear  dispersion  codes  whose  properties  are  intrinsically  governed  by  the 
Hurwitz-Radon  (HR)  families  of  matrices.  We  first  explore  the  quasi-orthogonal  codes  of  the  type  shown  in 
[11,25-28].  It  will  be  shown  that  all  4  x  4  quasi-orthogonal  codes,  where  each  entry  of  the  code  matrix  is  a 
symbol  variable  up  to  a  sign  change  and/or  complex  conjugation,  can  be  generated  from  any  two  independent 
codes  by  elementary  operations.  This  result  is  a  fundamental  unification  of  all  existing  (as  well  as  numerous 
previously  unrevealed)  4x4  quasi-orthogonal  codes  in  this  category.  If  all  symbols  are  selected  independently 
from  a  common  constellation  (which  will  also  be  referred  to  as  common  constellation  condition),  the  quasi- 
orthogonal  codes  may  have  diversity  two.  More  precisely,  half  of  the  quasi-orthogonal  codes  always  have 
diversity  two  under  the  common  constellation  condition,  and  the  other  half  may  have  diversity  either  two  or 
four  under  the  common  constellation  condition.  If  the  common  constellation  is  an  odd-numbered  phase-shift- 
keying,  half  of  the  quasi-orthogonal  codes  actually  have  diversity  four.  As  illustrated  later,  both  the  2x2 
orthogonal  codes  and  the  4x4  quasi-orthogonal  codes  can  be  expressed  as  linear  dispersion  codes  in  terms  of 
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the  HR  families  of  matrices.  This  observation  motivated  us  to  explore  8x4  linear  dispersion  codes  using  the 
HR  families  of  matrices  of  size  eight.  It  will  be  shown  that  a  class  of  8  x  4  linear  dispersion  codes  constructed 
with  the  HR  families  of  matrices  have  diversity  three  even  under  the  common  constellation  condition.  To  our 
knowledge,  this  is  the  first  8x4  unit-rate  linear  dispersion  code  that  is  guaranteed  to  have  diversity  three  when 
all  symbols  are  independently  selected  from  any  given  constellation.  This  is  an  unique  insight,  which  is 
unknown  from  the  previous  studies  of  linear  dispersion  codes  [19,22,29].  A  high  diversity  order  regardless  of 
symbol  constellation  is  practically  useful  since  it  could  reduce  the  physical  layer  complexity  associated  with 
constellation  constraint. 

In  Section  2,  we  review  the  HR  families  of  matrices.  In  Section  3,  we  show  that  all  4  x  4  quasi-orthogonal 
codes  can  be  constructed  by  two  independent  codes,  and  their  properties  are  discussed.  In  Section  4,  we 
introduce  a  class  of  8  x  4  linear  dispersion  codes  using  the  HR  matrices  of  size  eight,  and  show  that  these 
codes  have  diversity  three  under  any  given  constellation.  The  proof  of  the  diversity  three  property  is  a  lengthy 
part  of  this  paper.  We  hope  that  interested  readers  will  find  the  proof  theoretically  insightful  as  it  reveals 
detailed  structures  in  the  problem.  In  Section  5,  we  provide  a  simulation  example  to  illustrate  the  performance 
of  the  non-orthogonal  8x4  HR  code. 

1.1.  Notations  and  terminologies 

•  Lower  case  letters  are  used  for  scalars. 

•  Upper  case  letters  are  used  for  matrices. 

•  Underlined  lower  case  letters  are  used  for  vectors. 

•  In  normal  script,  *  denotes  an  undetermined  quantity.  As  superscript,  denotes  complex  conjugation. 

•  In  normal  script,  j  denotes  V— I.  In  subscript,  j  denotes  an  integer. 

•  Kronecker  product  is  denoted  by  0  as  defined  later. 

•  =  denotes  “equal  by  definition”. 

•  1 1  is  an  /  x  /  identity  matrix. 

•  As  superscripts,  T  denotes  transpose,  and  H  denotes  conjugate  transpose. 

•  91  denotes  real  part,  and  3  imaginary  part. 

•  Unless  specified  otherwise,  by  “symbol”,  we  mean  a  complex  variable. 

•  Unless  specified  otherwise,  two  vectors  of  variables  are  said  to  be  orthogonal  only  if  they  are  orthogonal  for 
all  values  of  the  variables. 

•  All  other  notations  are  defined  the  first  time  they  are  used. 

2.  HR  matrices 

2.1.  General  properties  of  HR  matrices 

Radon  Theorem  [30]:  Within  the  space  of  L  x  L  integer  matrices,  there  is  a  family  of  p  matrices 
{A0,Ai, . .  .,Ap_ i}  satisfying  A0  =  II  (the  L  x  L  identity  matrix)  and: 

•  Property  1(a):  AtAj  =  I L. 

•  Property  1(b):  At  =  —Aj  (/>  0). 

•  Property  2(a):  AjAj  =  —AjAt  (i^j), 

where  the  maximum  value  pmax  of  p  is  governed  by  L  as  follows.  Let  L  =  2 ab  where  b  is  odd,  a  =  4c  +  d  and 
0  ^  d  ^  3,  then  />max  =  8c  +  2d. 

A  family  of  matrices  defined  above  is  called  an  HR  family  of  matrices  of  size  L.  The  following  properties  of 
the  HR  matrices  follow  readily  from  Properties  1  and  2(a): 


•  Property  2(b):  AtAj  =  — AjAj  (i^j). 

•  Property  3:  For  any  real  vector  v,  vtAJAj  v  =  3(i  —  j) \v  \ 2  where  <5(x)  is  one  when  x  =  0  and  zero  otherwise, 
and  I  v  I  is  the  norm  of  v. 
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•  Property  4:  Given  the  matrices  Aj,j  =  1, ...  ,7V,  within  an  HR  family,  we  have 

(N  \  (  N  \  T  N 

It  XJAJJ  aJAJJ  =  U  XP ■ 


An  HR  family  of  matrices  of  any  size  can  be  constructed  from  the  following  2x2  matrices  [30]: 

0  1 


P=  [  °  1 V  Q=(l  0 
loy  ^  Vo  -l 


R  = 


-1  0 


(3) 


2.2.  HR  matrices  of  size  two 

For  L  =  2,  an  HR  family  is  {I2,R}.  The  Alamouti  code  can  be  expressed  in  terms  of  the  2x2  HR 
matrices,  e.g., 


C(si,s2)  = 


/x2(l)\  fxflf 
vi(2)  )  +  x2(2)  )’\  *i(2) 


+  j* 


'*2(1)' 

*2(2) 


/  xi(2)  +  jx2(l)  vi(l)+jx2(2)\ 

\  — *t(l)  +  j*2(2)  vi(2)-jx2(l)  y 
where  s\  =  xi(2)  +  jx2(l)  and  s2  =  vi(l)  +  jx2(2). 

2.3.  HR  matrices  of  size  four 

For  L  =  4,  an  HR  family  consists  of  the  following  matrices: 
Qo  =h,  QX=P®R,  Q2  =  R®I2 ,  Q3  =  Q®R, 
where  ®  is  the  Kronecker  product,  e.g., 

fan  aX2  \  B  _  ( anB 
l  a2\  a22  I  \  a2\B  a22B  I 


(4) 


(5) 


The  HR  families  of  size  four  are  closely  related  to  4  x  4  quasi-orthogonal  codes  as  discussed  later. 

The  following  theorem  provides  the  complete  set  of  HR  families  of  matrices  of  size  four.  If  the  first  matrix  in 
each  family  is  fixed  to  be  identity,  there  are  total  2  x  23  =  16  HR  families  of  size  four. 

Theorem  2.1.  Any  HR  family  of  matrices  of  size  four  has  either  one  of  the  following  two  possible  forms : 

fii={Go;±&;±fi2;±03}  and  ^2={g0;±gi;±g2,±g3}, 

where  Go  =  Q0,  G\  =  QX[Q  0  Q],  G2  =  Q2[—I2  0  Q],  and  G3  =  Q3(Q  ® /2).  The  Qt  matrices  were  defined 
previously.  {Note  that  in  this  paper ,  ±  in  one  place  should  be  treated  as  independent  of  ±  in  another  place  unless 
specified  otherwise .) 

Proof.  The  proof  of  this  theorem  requires  an  exhaustive  but  finite  search,  which  is  tedious  but  feasible.  In  the 
following,  we  provide  an  outline  of  the  proof.  The  goal  is  to  show  that  under  Properties  1(a),  1(b)  and  2,  only 
Q i  and  Q2  can  be  valid  HR  families. 

Let  us  first  search  for  all  possible  4x4  HR  matrices  satisfying  Properties  1(a)  and  1(b).  Under  Property  1(a), 
each  entry  of  an  HR  matrix  F  is  either  zero  or  ±  1  and  each  row  of  F  has  no  more  than  one  non-zero  entry. 
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Under  Property  1(b),  i.e.,  FT  =  —F  where  FVU,  we  can  write 

/  B\,\  B\  2  \ 
yi?2,i  B22  j 

where  Bq  is  a  2  x  2  block  matrix.  It  follows  that  Bji  =  —Bij  for  i  =  1,2  and  Bj2  =  —B2\.  Therefore,  B\  \  and 
B22  are  both  equal  to  either  the  2x2  zero  matrix  or  ±R,  where  B\  \  and  B2  2  may  have  the  same  or  opposite 
signs.  Namely,  we  have 


either 


B\,\  =  B22  =  0, 
Bh  =  ~BV  *0 


or 


f  B i,i  =  ±R, 

J  52, 2  =  ±R, 

B\,2  =  B2  \  =  0. 


(6) 


It  is  easy  to  verify  that  if  B\j  =  B22  =  0,  the  only  possible  choices  of  B\  2  are  ±/2,±P,  ±R,  ±Q. 

By  an  exhaustive,  finite  and  feasible  search,  one  can  further  verify  that  any  HR  matrix  F  of  size  four  has  to 
be  one  of  the  matrices  in  Q\  or  Q2.  In  other  words,  the  above  two  sets  Qi  and  Q2  contain  all  possible  4x4  HR 
matrices  satisfying  Properties  1(a)  and  1(b). 

We  next  consider  Property  2(a).  It  is  straightforward  to  verify  that  Q \  and  Q2  each  satisfies  Property  2(a). 
Namely,  each  of  Q  \  and  Q2  (with  a  fixed  set  of  ±  signs)  is  a  valid  HR  family.  In  order  to  prove  that  there  is  no 
other  possible  HR  family,  we  need  to  show  that  for  any  i,j  with  l<ij <3,  Qt  and  Gj  cannot  co-exist  in  one 
family.  Indeed,  it  is  straightforward  to  verify  that  QjGj  ^  —  G]Qi,  be.,  Property  2(a)  is  not  satisfied  by  Qt  and 
Gj  for  any  i,j  with  1  ^z,y'^3.  □ 


2.4.  HR  matrices  of  size  eight 

For  L  =  8,  an  HR  family  is  determined  by  the  following  eight  matrices  [30]: 

hJ2®  R®  I2,I2®  P  ®  R\Q®  Q®R\P®  Q®R\R®P®  Q]R®P®P\R®  Q®I2  (7) 

which  is  easy  to  verify. 

An  HR  family  of  matrices  of  size  eight  (or  size  integer-power- of- two  no  less  than  eight)  have  the  following 
Properties  5  and  6: 

•  Property  5:  Given  distinct  i,  k,  m,  n,  then  AjAmAjAn  is  unitary  and  symmetric.  And  the  eigenvalues  of 
AjAmAjAn  are  ±1  of  differing  signs. 

Proof.  It  is  easy  to  verify  that  AjAmAjAn  is  unitary.  To  prove  the  symmetry,  we  apply  Property  2: 

Aj  A m A  ■  A n  =  AmAjAnAi 

— 

=  —  AnAmAj  Aj 
=  AlAtAmAi 

=  (A]AmAjA„)T.  (8) 

We  now  prove  the  eigenvalue  property.  The  symmetry  and  orthogonality  of  AjAmAjAn  tell  us  that  each  of  its 
eigenvalues  is  +1  or  —1.  We  only  need  to  prove  that  the  signs  of  these  eigenvalues  are  always  mixed.  Suppose 
that  the  eigenvalues  of  AjAmAjAn  are  all  equal  to  one  or  all  equal  to  minus  one.  Then,  the  matrix  is  either  /  or 
-I,  i.e., 

AjAmAjAn  =  ±IL 
=>  AjAn  =  ±AjnAj 

=>  An  =  ±AiATmAj,  (9) 
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where  the  signs  in  the  above  three  equations  are  consistent  with  each  other.  Considering  another  member  Ai  in 
the  HR  family  where  /  is  distinct  from  m,  n,  ij ,  we  have 

A„Aj=  ±  AjA^AjAj 

=  =f A,AjAmAj 

=  =F A,ATmAjAj 

=  ±A,A]AmAj,  (10) 

where  the  signs  in  the  above  equations  are  consistent  with  each  other.  From  (9),  we  have 

-AiAl  =  =p AiAjAmAj .  (1 1) 

Since  AnAj  =  —AiAjt,  (10)  and  (11)  imply  —AiAjAmAj  =  AiAjAmAj  and  hence  A/AjAmAj  =  0.  This 
contradicts  the  condition  of  non-zero  eigenvalues.  Therefore,  all  eigenvalues  of  AjAmAjAn  are  ±1  of  differing 
signs.  □ 

•  Property  6:  If  AjAkAjnAn  ±  AjAtA^Ar  =  0  where  all  matrices  are  distinct  HR  matrices  in  one  family, 
we  have 

AlAk,AlAn,  ±  (~\)E A]AhA]An  =  0, 

where  E  is  an  integer  and 

[Ai\ ,  Ak  j ,  Ami ,  Ani ,  Ajx ,  Atl ,  Aqx ,  Ari  ] 

is  a  matrix  series  generated  by  exchanging  E  pairs  of  matrices  in  the  following  matrix  series: 

[At,  Ak,  Am,  An,  Aj,  At,  Aq,  Ar\. 


Proof.  It  suffices  to  prove  that  by  switching  any  pair  of  matrices  in  AjAkAjnAn  ±  AjAtA^Ar  =  0,  the  ±  sign 
changes  to  the  qp  sign. 

First,  we  consider  the  two  matrices  exchanged  are  from  either  the  first  term  AjAkAjnAn  or  the  second  term 
AjAtA^Ar.  Exchanging  any  two  matrices  A  and  B  that  have  d  other  matrices  in  between  is  equivalent  to 
the  following  process:  (a)  exchanging  A  with  its  next  matrix  repeatedly  until  A  is  right  behind  B  and  then  (b) 
exchanging  B  with  its  preceding  matrix  repeatedly  until  B  is  in  the  original  position  of  A.  Step  (a)  undergoes 
d  +  1  exchanges  of  neighboring  matrices,  and  step  (b)  undergoes  d  exchanges  of  neighboring  matrices.  By 
Property  2(a),  the  resulting  expression  has  changed  its  sign  2d  +  1  times,  and  therefore  there  is  a  net  sign 
change. 

We  now  need  to  prove  the  Property  6  for  the  case  where  a  pair  of  matrices  exchanged  are  “crossed”  between 
the  two  terms.  We  first  observe  the  following  equivalent  expressions: 


AjAkATmAn±A]AtATqAr  =  0, 

AjAiAjAkAjnAn  ±  AjAiAjAtAJqAr  =  0 
AjAkAjnAn  =f  AjAtA^Ar  =  0, 


(12) 

(13) 


where  we  should  notice  that  only  the  leading  matrices  A{  and  Aj  are  actually  exchanged  when  we  move  from 
(12)  to  (13).  We  now  consider  a  matrix  A  (from  the  first  term)  that  has  dA  matrices  proceeding  it  and  another 
matrix  B  (from  the  second  term)  that  has  ds  matrices  proceeding  it.  In  order  to  exchange  A  and  B,  we  can 
do  the  following  steps:  (a)  move  A  to  the  front  of  the  first  term  and  move  B  to  the  front  of  the  second  term, 
(b)  exchange  A  and  B,  and  (c)  move  B  to  where  A  initially  was  and  move  A  to  where  B  initially  was.  We  see 
that  step  (a)  undergoes  the  sign  change  dA  +  dB  times,  step  (b)  undergoes  the  sign  change  once,  and  step  (c) 
undergoes  the  sign  change  dA  +  dB  times.  Therefore,  there  is  a  net  sign  change.  □ 
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Given  the  HR  family  of  size  eight  shown  in  (7),  it  can  be  verified  that 

AjA6AjA0  =  AjA2AjA3.  (14) 

(At  this  stage,  it  is  not  clear  whether  (14)  holds  under  a  more  general  condition.)  Let  i,k,m,n,j,t,q,r  be  a 
permutation  of  [0, 1, ,  7].  Then,  together  with  Properties  6  and  3,  (14)  implies  the  following  properties: 


•  Property  7(a):  A[AiAjnAn  =  ±AjAtA^Ar  where  the  sign  depends  on  the  ordering  of  the  indices. 

•  Property  7(b):  A^A„  =  =fAjAfAjAtAjAr  where  the  sign  is  consistent  with  Property  7(a). 

•  Property  7(c):  For  any  real  vector  v,  vT(AjA,-AjAtAjAr)v  =  0. 

3.  Quasi-orthogonal  codes 


We  now  present  a  complete  family  of  4  x  4  quasi-orthogonal  codes  (or  code  matrices)  of  Type  I.  A  4  x  4 
quasi-orthogonal  code  matrix  of  Type  I  is  defined  as  such  that  each  entry  in  the  matrix  is  an  element  from 
the  symbol  set  {±s^,±s^f±s^,±s^}  and  a  pair  of  columns  of  the  matrix  is  orthogonal  to  the  other  pair 
(and  the  two  columns  in  each  of  the  above  pairs  are  not  necessarily  orthogonal  to  each  other).  Note  that  each 
±  is  an  independent  plus  or  minus  sign,  and  each  superscript ( denotes  independently  the  presence  or  absence 
of  complex  conjugation.  The  above  definition  of  quasi-orthogonal  code  of  Type  I  was  used  in  [11].  It  is 
obvious  that  if  £(51,52,53,54)  is  a  Type  14x4  quasi-orthogonal  matrix  of  the  four  symbols  51,52,53,54,  then 
numerous  Type  I  quasi-orthogonal  codes  can  be  constructed  by  the  following  elementary  operations: 


00(1,52,^4)  =  P'S(±$,±$MgMg)Pc, 


(15) 


where  (k\,k2,k'i,kf)  is  a  permutation  of  (1,2, 3,4),  PT  permutes  the  rows  and/or  reverses  the  signs  of  none  or 
some  rows,  and  Pc  permutes  the  columns  and/or  reverses  the  signs  of  none  or  some  columns.  Note  that  the 
above  statement  is  obvious  because  none  of  the  operations  Pr,  Pc,  ±  and  *  changes  the  quasi-orthogonality  of 
(15).  While  the  above  statement  is  obvious,  a  number  of  Type  I  quasi-orthogonal  code  matrices  have  been 
introduced  in  the  literature  without  mentioning  the  connections  among  them.  For  beginners,  each  of  these 
codes  appears  to  be  a  new  one.  Even  for  experts,  it  was  unknown  whether  all  existing  codes  of  Type  I  are 
related  to  each  other  by  (15).  In  this  section,  we  will  show  that  not  all  Type  I  codes  are  related  to  each  other  by 
(15),  but,  however,  there  are  only  two  groups  of  Type  I  codes.  The  codes  in  each  of  the  two  groups  are  related 
to  each  other  by  (15),  but  no  code  from  one  group  is  related  to  any  code  from  the  other  group. 

One  can  also  extend  the  family  of  quasi-orthogonal  codes  by  allowing  left  or  right  multiplication  of  a 
diagonal  matrix  to  (15),  which  is  a  simple  extension  of  the  Type  I  codes.  There  are  also  orthogonal  or  quasi- 
orthogonal  codes  where  the  entries  of  the  code  matrix  are  non-linear  functions  of  the  symbol  vector 
{51,52,53,54}  and/or  the  symbol  constellation  is  constrained  [31].  In  this  paper,  we  will  not  consider  any  quasi- 
orthogonal  codes  other  than  the  Type  I  quasi-orthogonal  codes. 

Two  codes  will  be  called  independent  of  each  other  if  they  are  not  related  to  each  other  according  to  (15),  or 
otherwise  dependent  on  each  other.  It  is  clear  that  a  complete  set  of  quasi-orthogonal  codes  can  be  generated 
by  all  independent  codes  via  (15).  But  we  will  show  that  the  number  of  independent  codes  is  two.  Examples  of 
such  two  independent  codes  are  also  given.  All  existing  codes  of  this  type  will  be  explicitly  expressed  in  terms 
of  the  two  independent  codes. 


3.1.  Independent  quasi-orthogonal  codes 

We  show  next  that  there  are  only  two  independent  4x4  Type  I  quasi-orthogonal  codes.  But  first,  we  have 
the  following  property  about  complex  vectors: 

Lemma  3.1.  Let  s  —  sr  +  j5f  be  a  4  x  1  complex  vector  in  the  four-dimensional  complex  space  C4  where  sr  is  the 
real  part  and  st  is  the  imaginary  part.  Define  the  second  4x1  complex  vector  as  p  =  Mrsr  +  j  where  Mr  and 
Mi  are  unitary  integer  matrices.  Then ,  sH  p  =  0  holds  for  all  s  in  C4  if  and  only  if  (to  be  denoted  by  iff)  a  pair  of 
elements  in  s  is  orthogonal  to  the  corresponding  pair  in  pfor  all  s  in  C4  and  the  other  pair  of  elements  in  s  is 
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orthogonal  to  the  other  corresponding  pair  in  p  for  all  s  in  C4.  (The  pair-wise  orthogonality  is  equivalent  to  that  in 
Alamouti  code  [1].) 

Proof.  Taking  the  real  and  imaginary  parts  of  =  0  separately,  we  have  =  0  iff  t\  +  ti  =  0  and  ti  =  0 
where  t\  =  sj:Mrsr,  t2  =  sj  and  t2  =  sf  Mrsr  —  sf  Mjsr. 

Because  of  the  independence  between  sr  and  sp  t\  +  t2  =  0  iff  t\  =  0  and  t2  =  0.  Because  both  sr  and  st  are 
any  real  vectors,  we  have  t\  =  0  iff  Mr  =  —  Mj  and  t2  =0  iff  Mt  —  —Mj.  Similarly,  we  have  T  =  0  iff 
Mr  =  Mj. 

The  above  implies  that  sHp  =  0  holds  for  all  s  in  C4  iff  Mt  =  —Mr  and  Mr  =  —Mj. 

With  the  above  choice  of  Mr  and  Mif  =  0  is  equivalent  to  sf  =  0.  We  next  apply  that  Mr  =  —Mj 
and  each  row  of  Mr  has  only  one  non-zero  entry  ±1.  Without  loss  of  generality,  we  can  assume  that  the  (z'o,y0) 
entry  of  Mr  is  ±1,  i.e.,  Mr(/0j0 )  =  ±1.  Then  Mr(j0,io )  =  -Mr(/0,y0)  where  io^j0  and  all  other  elements  in  the 
/oth  and y0th  rows  of  Mr  are  zero.  This  property  implies  that  the  /oth  and y0th  elements  in  s  cancel  each  other  in 
the  form  ±(s*(io)s*(j0)  —  s*(jo)s*(io))in  sfMrs*  =  0.  In  other  words,  the  /oth  and  y0th  elements  in  s  are 
orthogonal  to  the  /oth  and  y0th  elements  in  p  for  all  s  in  C4.  More  explicitly,  the  /oth  and y0th  elements  in  p  are 
±?(jo)  and 

Following  the  same  reasoning,  p  =  0  for  all  s  in  C4  iff  the  other  two  elements  in  s  are  also  orthogonal 
(in  the  same  way  as  described  above)  to  the  other  two  corresponding  elements  in  p  for  all  s  in  C4.  □ 

From  Lemma  3.1,  the  next  theorem  follows  (which  corrects  a  result  shown  in  [32]): 

Theorem  3.1.  Define  a  code  set  Sq  of  4  x  4  quasi-ortho gonal  codes  for  the  symbol  vector  s  =  (s\,S29s$,sfF , 
where  each  element  in  a  (normally  full  rank )  code  matrix  has  the  form  ±s^  and  two  of  the  four  columns  in  each 
code  are  orthogonal  to  the  other  two  over  all  s  in  C4.  Then,  the  following  two  codes  S\(s\,s2,s$,sf)  and 
S2(Si,S2,S3,S4): 


(  S\ 

-S4 

r.* 

S1 

-4\ 

(  Si 

s4 

4 

-4\ 

Si 

S  3 

_  0* 

S\ 

-4 

Si 

Si 

-4 

S4 

Si 

-si 

_c* 

S4 

Sl 

,  S2  = 

Sl 

Si 

_c* 

S4 

4 

,S4 

S\ 

c* 

£3 

S1  ) 

iS4 

Si 

S1 

—Sj  y 

span  all  the  codes  in  S q  via  (15),  and  furthermore  S\{s\9s2,s^,sf)  and  S2(s\9S2,si9sf)  are  not  related  to  each 
other  via  (15). 

Proof.  We  will  say  that  a  pair  of  codes  are  dependent  of  each  other  if  they  are  related  via  (15),  or  otherwise 
independent  of  each  other.  Our  proof  is  constructive  in  that  we  will  construct  a  largest  possible  set  S  of 
independent  codes.  It  is  important  to  stress  that  permutations  of  rows  and/or  columns,  permutations  of  symbol 
indices,  change  of  sign  to  each  row  and/or  column,  and  sign  and/or  conjugation  changes  to  each  symbol  are  all 
variations  allowed  by  (15)  among  dependent  codes.  Our  proof  consists  of  several  steps  by  which  the  above 
variations  are  eliminated  from  a  largest  possible  set  of  independent  codes.  These  steps  will  lead  to  two  possibly 
independent  codes  S\  and  S2.  These  two  codes  Si  and  S2  are  then  finally  verified  to  be  independent. 

Without  loss  of  generality,  we  can  fix  the  first  column  of  each  code  matrix  in  S  to  be 
Furthermore,  we  can  choose  S  in  such  a  way  that  the  first  two  columns  of  each  code  matrix  are  orthogonal  to 
the  last  two  columns.  From  Lemma  3.1,  it  follows  that  among  all  code  matrices  in  S,  there  are  no  more  than 
the  following  two  possible  forms  up  to  the  variations  defined  by  (15): 


(  Si 

* 

c* 

S1 

*\ 

(  Si 

* 

-S*2 

Si 

* 

-s*i 

* 

Si 

* 

s\  * 

or  T2  = 

Si 

* 

-sl 

* 

Si 

* 

s4  * 

yS4 

* 

0* 

*V 

^  s4 

* 

^3  *y 

where  *  (not  in  superscript)  denotes  a  unspecified  entry.  Note  that  in  each  of  T \  and  T2,  the  first  two  elements 
of  the  first  column  are  orthogonal  to  the  first  two  elements  of  the  third  column  for  all  s  in  C4,  and  the  last  two 
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elements  of  first  column  are  orthogonal  to  the  last  two  elements  of  the  third  column  for  all  s  in  C4.  One  can 
easily  verify  that  no  other  possibility  exists  that  is  independent  of  T\  and  T2. 

Similarly,  from  T\,  one  can  generate  no  more  than  the  following  four  possibilities  in  S  up  to  (15): 


(  Si 

* 

s* 

-4) 

(  Si 

* 

S2 

-4) 

S2 

* 

-4 

_c* 

S4 

,  Tia  = 

S2 

* 

-4 

* 

s4 

S3 

* 

-sj 

Sf 

S3 

* 

-S4 

Si 

^4 

* 

st 

s|  y 

^s4 

* 

o* 

S3 

-s\ 

(  Si 

* 

sj 

-sj\ 

(  Si 

* 

S2 

-sj\ 

S2 

* 

—  V* 

Si 

jk 

—  c* 

S3 

sk 

,  T 154  = 

S2 

* 

—  V* 

Si 

>k 

S3 

S3 

* 

— sj 

s2 

S3 

* 

— sj 

— sX 

^  S4 

* 

0* 

S3 

C** 

/ 

^  S4 

* 

0* 

S3 

Si  J 

Furthermore,  from  T 2,  one  can  generate  no  more  than  another  four  possibilities  in  §  up  to  (15): 


(  Si 

* 

_c* 

S2 

-4\ 

(  Si 

* 

_c* 

S2 

-4\ 

S2 

* 

Si 

4 

S2 

* 

c* 

Sl 

—  c* 

S4 

S3 

* 

_c* 

s4 

4 

,  T  2,2  — 

S3 

* 

_o* 

s4 

4 

^  S4 

* 

S3 

~s2  J 

^  S4 

* 

0* 

S3 

4  ) 

(  Si 

* 

~S2 

-4) 

(  Si 

* 

~S2 

-4) 

S2 

* 

st 

-4 

S2 

* 

st 

4 

* 

„ik 

,  T2, 4  =3 

* 

„ik 

S3 

-sj 

S2 

S3 

“Sj 

—s\ 

^s4 

* 

o* 

S3 

Sl  j 

^s4 

* 

o* 

S3 

Sl  j 

By  filling  the  second  column  of  each  of  the  above  matrices  (to  satisfy  the  orthogonality  condition),  it  follows 
that  there  are  no  more  than  the  following  eight  possibilities  in  S  up  to  (15): 


(  Sl 

S4 

s\ 

-4\ 

(  Si 

s4 

s\ 

-s?\ 

S2 

S3 

-St 

_o* 

S4 

S2 

S3 

-st 

Sj 

Tu  = 

S3 

-S2 

-S4 

4 

,  Tu  2  = 

S3 

S2 

_  o* 

s  4 

st 

’ 

y  S4 

Sl 

St 

4  ) 

y  S4 

Sl 

c* 

S3 

~S*2  } 

(  Sl 

S3 

s\ 

-4) 

(  Si 

S3 

S2 

— sj 

:\ 

S2 

-S4 

-St 

-4 

S2 

s4 

-st 

st 

Tl,3  = 

* 

,  Th4  = 

S3 

Sl 

-S4 

4 

S3 

-Sl 

-S4 

-s\ 

^S4 

s2 

0* 

S3 

4  ) 

^s4 

— S2 

o* 

S3 

o* 

Sl 

/ 

(  Sl 

s4 

_o* 

S2 

-4\ 

(  Sl 

-S4 

_ 0* 

S2 

_  0* 

S3 

;\ 

S2 

S3 

Sl 

S4 

S2 

S3 

0* 

Sl 

— S4 

T  2,1  = 

,  t2, 2  = 

jk 

S3 

-S2 

^4 

st 

S3 

S2 

“Sj 

st 

^  S4 

-Sl 

S3 

-S|  y 

^  S4 

-si 

S3 

sl 

/ 
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(  S\ 

^3 

-s\ 

-sj\ 

(  Si 

S3 

-S2 

-Sj\ 

S2 

-S4 

s\ 

-4 

S2 

s4 

S* 

^3 

S3 

~S\ 

_o* 

S4 

O* 

^2 

,  T2A  = 

S3 

Si 

~S*4 

— S2 

^  S4 

S2 

S3 

Si  j 

^  S4 

s2 

S3 

Si  j 

If  we  let  *Si  =  T\t\  and  S2  =  7\2,  then  it  can  be  verified  that 


Thi(s)  =  P\\S2(-si,s2,S4,S3)P2\,  (17) 

TiM  =  PnSi(s3,S2,S4,si)P22,  (18) 

T2,i(sJ  =  PnSi(s4,  s2,  S3,  s\)P23,  (19) 

T2,2(£)  =  P\4S2{-S\,S2,S3,S4)P24,  (20) 

P  2,3  (s)  —  Pl5Sl(sUS2,S4,S3)P25,  (21) 

T  2,4(1)  —  Pi6S2(si,  S2,  S4,  S3)P26,  (22) 

where 

Tfil  =  ^-(1,1), +(2,2), +(3,4), +(4,3),  (23) 

Pl\  =  P +(1,1), — (2,2), — (3,3), — (4,4)  5  (24) 

P\2  =  ^+(1,4), +(2,2), +(3,1), +(4,3),  (25) 

P22  =  7?(l,l),+(2,2),-(3,4),-(4,3),  (26) 

P 13  =  7J+(l,4),+(2,2),+(3,3),+(4,l),  (27) 

7^23  =  75+(l)l))+(2,2),— (3,4),— (4,3) ,  (28) 

7>14  =  ^-(1,1), +(2,2), +(4,3), +(3,4),  (29) 

7^24  =  75+(i;i))+(2,2),+(3,3),— (4,4),  (30) 

7*15  =  7>+(l,l),+(2, 2), +(4,3), +(3,4),  (31) 

P25  =  7>+(l,l)j— (2,2),— (3,3),+(4,4),  (32) 

P 16  =  7?+(i)i))+(2,2),+(4,3),+(3,4),  (33) 

P2  6  =  7>+(l)l),+(2,2),-(3,3),+(4,4)  (34) 


and  7?±0l!7!)!±(/2s/2)5±(/3!/-3X±(/4!/4)  denotes  a  matrix  where  the  entries  at  0'i  Ji),  O2J2),  (*3j3),  O4J4)  are  ±1  and  all 
other  entries  are  zero. 

Therefore,  there  are  no  more  than  two  independent  codes  in  S,  and  -Si  and  S2  are  two  possible  independent 
codes. 

To  prove  that  Si  and  S2  are  indeed  independent  codes,  we  need  to  observe  a  property  from  (15).  If  all 
elements  of  s  come  from  a  common  symmetric  constellation  (symmetric  in  terms  of  sign  change  and  complex 
conjugation),  then  it  is  obvious  from  (15)  that  any  two  dependent  codes  Ci(s)  and  C2(s)  must  satisfy  the 
following  identity: 

minrankfCifs;))  =  minrank(C2(|)),  (35) 

5+0  l+o 
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where  s  is  a  function  of  s,  satisfying 

i=(±£\±$,±$M£)T. 

In  the  next  subsection,  we  show  that  *Si  and  S2  do  not  satisfy  the  property  (35).  In  fact,  we  will  consider  an 
equivalent  situation  where  s  is  replaced  by  A  s  and  the  constellation  of  A  s  is  symmetric  with  respect  to  sign 
change  and  complex  conjugation.  Since  S 1  and  S2  do  not  satisfy  the  necessary  condition  (35)  as  required  for 
any  pair  of  dependent  codes,  *Si  and  S2  are  independent.  □ 


3.2.  Diversity  of  quasi-orthogonal  codes 


The  diversity  of  a  code  matrix  C(s)  is  the  minimum  rank  of  C(A  s)  over  all  possible  A  s  9^  0  where  A  s  is  the 
difference  between  two  symbol  vectors.  For  the  diversity  analysis  of  quasi-orthogonal  codes,  we  will  assume 
that  the  constellation  of  A  s  is  symmetric  with  respect  to  sign  and  conjugation.  Such  a  condition  is  common  in 
practice.  Because  of  the  property  (35)  among  dependent  codes,  to  study  the  diversity  of  all  the  4  x  4  quasi- 
orthogonal  codes,  it  suffices  to  consider  the  diversity  of  Sfs)  and  S2(s). 

The  diversity  of  Sfs)  is  the  minimum  rank  of  A-(Ay)  over  all  A s  ^0  where  As  =  (A^i,  Aj2,  AJ3,  A^)1.  Here, 


Di(s£=Si(s)aSl(s)  = 


and 


D2(sJ=S2(s)hS2(s)  = 


where 


a(s)  ]bi(s)  0  0 

— j^iGO  a(f)  0  0 

0  0  a(s)  }bi(s) 

0  0  -j  bi(s)  a(s) 

a(s )  b2(s)  0  0 

b2(s)  a(s )  0  0 

0  0  a(s)  b2(s)  ’ 

0  0  b2(s)  a(s) 


(36) 


(37) 


(38) 


bi(s)  =  23(sis*4  +  5^3),  (39) 

b2(s)  =  2SR(jIj4  +  ^3).  (40) 

It  is  clear  from  (36)  and  (37)  that  the  rank  of  D/(As)  is  either  two  or  four  as  long  as  A s  ^0.  To  determine 
the  conditions  under  which  A(A;S)  has  the  full  rank  (i.e.,  rank  four),  it  is  useful  to  examine  its  determinant. 
The  determinants  of  Dj{s )  are  given  by  det(A(y))  =  («(^)  -  bi(s))2(a(s)  +  bi(s))2  or  equivalently 

det(D\(s))  =  [|jji  -54|2  +  \s2  -  j^l2]2[lj^i  +^4l2  +  1^2  +j^3l2]2,  (41) 

det(Z)2(j))  =  [\si  +  54|2  +  \s2  +  s3\2]2[\si  -  ^4|2  +  \s2  -  s3\2]2.  (42) 


We  show  next  that  the  conditions  for  a  full  rank  A  (A  s)  is  generally  different  from  that  for  a  full  rank  D2( A  s) 
even  if  the  constellation  of  A  s  is  symmetric.  (Note  that  this  was  needed  to  complete  the  proof  of  Theorem  3.1.) 

Expression  (41)  would  make  Lemma  2.1  in  [33]  more  complete  and  would  also  enrich  Lemma  2.2  in  the 
same  paper.  It  follows  from  (41)  that  if  and  only  if  the  set  {jA^i}  and  the  set  {±A^4}  have  no  common  elements 
except  zero  or  the  set  {A^}  and  the  set  {ijA^}  have  no  common  elements  except  zero,  then  A(A5)  has  full 
rank  as  long  as  As  ^0. 

Expression  (42)  is  given  and  well  discussed  in  Lemma  2.1  in  [33].  It  follows  from  (42)  that  if  and  only  if  the 
set  { A^i }  and  the  set  {±As4}  have  no  common  elements  except  zero  or  the  set  {Ay?}  and  the  set  {±A^3}  have  no 
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Fig.  1.  The  set  of  A st  (crosses)  and  the  set  of  jA.s;-  (circles)  for  3-PSK.  Note  that  the  constellation  of  As,  is  symmetric  in  terms  of  sign  change 
and  complex  conjugation,  which  meets  the  condition  for  (35). 

common  elements  except  zero,  then  D2(A  s)  has  full  rank  as  long  as  A  s  ^  0.  This  is  an  observation  also  made  in 
[13-16]. 

For  example,  if  the  constellation  is  odd-numbered  phase-shift-keying  as  illustrated  in  Fig.  1,  then  the  set 
{jAsi}  and  the  set  {iA^}  have  no  common  elements  except  zero  and  the  set  {A^}  and  the  set  {ijA^}  have  no 
common  elements  except  zero,  and  hence  Z>i(A5)  has  full  rank  under  As  ^0.  But  for  the  same  constellation, 
D2(As)  does  not  have  full  rank  under  As  ^0. 

3.3.  Quasi- orthogonal  codes  in  terms  of  HR  matrices 

We  now  illustrate  that  all  quasi-orthogonal  codes  are  linear  dispersion  codes  by  expressing  the  two 
independent  codes  S\  and  S2  in  terms  of  the  HR  matrices.  Let  the  real  and  imaginary  parts  of  each  symbol  Sk 
be  expressed  as  Sk  =  rjc  +  j h-  It  is  not  difficult  to  verify  the  following  results.  For  the  first  code, 


(r\ 

-r4 

-r2 

-^3\ 

(h 

—i4 

h 

h  \ 

r3 

r\ 

-r4 

h 

h 

-h 

U 

Si  = 

r3 

—r2 

r4 

r\ 

+  j 

h 

-h 

~i4 

~i\ 

\r4 

r\ 

~r- i 

r2  ) 

V'4 

i\ 

h 

-h) 

—  [QqLx,  —Q\L\,  Q2H.\ •>  —Q2L\\  +  ][Qih->  -Q3IU  ~Qih>  Qoh\> 

where  r_x  =  [rx, . . .  ,r4]T  and  ix  =  [z'3,/4,  —i\,  —i2]J.  For  the  second  code, 


/  r\ 

r4 

^2 

-r2\ 

/  h 

i4 

—  h 

h  \ 

r2 

?3 

-r\ 

r4 

h 

h 

h 

—i4 

s2  — 

r3 

r4 

r\ 

+  j 

h 

h 

i4 

~i\ 

U 

n 

r3 

-r2) 

u 

h 

-h 

h  / 

—  [—Q3KL2,  Q2r2,  Q0Kr2,  Qjr2 ]  +  j[—  QqKi2,  Qji2,  Q2Ki2,  Q2i2],  (43) 

where  r2  =  \-r2,  -ri,r4,r3]T,  i2  =  [iu  -i2,i3,  -i4]T,  and  K  =  —I2  0  Q. 
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Table  1 

Expressions  of  existing  quasi-orthogonal  codes  in  terms  of  two  independent  quasi-orthogonal  codes 


Authors/reference  Original  code  PT  Si  Pc 


Papadias  and  Foschini  [25] 

(  ^ 

Sl 

S3 

s4 

\ 

+(h3) 

+(2,1) 

(4\ 

+(U) 

-(2,4) 

4 

-4 

V* 

s4 

~4 

+(4,4) 

+(3,2) 

Si 

S3 

-(3,2) 

-(4,3) 

S3 

-S4 

-Sl 

Si 

Sl 

\4 

4 

-4 

-4 

/ 

{4) 

Hou  et  al.  (20)  in  [26] 

(  Sl 

Sl 

S3 

S4 

+(1,1) 

+(2,4) 

(  SX 

\  +(U) 

+(2,2) 

-si 

Sl 

-S4 

S3 

+(3,3) 

+(4,2) 

Si 

-(3,4) 

+(4,3) 

-4 

4 

4 

-4 

-*3 

[-4 

-4 

4 

4 

) 

\-S2, 

/ 

Ran  et  al.  (10)  in  [27] 

/Si 

Sl 

S3 

s4\ 

+(U) 

+(2,2) 

/si\ 

+(U) 

+(2,4) 

4 

-4 

-4 

S3 

+(3,3) 

+(4,4) 

Si 

si 

+(3,2) 

-(4,3) 

4 

-4 

-4 

S3 

V  s^ 

S3 

Sl 

Sl/ 

\S4/ 

Tirkkonen  et  al.  (10)  in  [28] 

(  si 

Sl 

S3 

s4\ 

+(U) 

+(2,2) 

/  Sl  ' 

\  +(U) 

-(2,3) 

4 

~S*4 

S3 

+(3,4) 

+(4,3) 

Si 

-si 

-(3,2) 

+(4,4) 

S3 

S4 

Sl 

S2 

-S4 

(-j 

V* 

S3 

~4 

4) 

V  S3  , 

/ 

Jafarkhani  [11] 

11 

Sl 

S3 

S4' 

\ 

-(1,1) 

+(2,2) 

/-Sl' 

\  +(U) 

-(2,4) 

~S2 

4 

~4 

4 

+(3,3) 

+(4,4) 

Sl 

-4 

+(3,2) 

-(4,3) 

~4 

~4 

Sl 

4 

-4 

\  S4 

S3 

— S2 

Si . 

/ 

\  S4  , 

/ 

A  column  vector  of  symbols  is  used  instead  of  a  row  vector  for  a  more  compact  form.  The  arrays  of  ±(zj)  in  the  third  and  fifth  columns 
indicate  the  ±1  entries  in  the  integer  matrices  Pt  and  Pc.  Other  unspecified  entries  of  P,  and  Pc  are  zero. 


3.4.  Previously  published  quasi-orthogonal  codes 

The  previously  published  4x4  unit-rate  quasi-orthogonal  codes  of  Type  I  can  all  be  expressed  in  terms  of 
the  above  two  independent  codes  via  (15).  Table  1  summarizes  these  connections. 


4.  Non-orthogonal  HR  codes 

In  this  section,  we  present  an  8  x  4  unit-rate  code  and  prove  that  it  has  diversity  no  less  than  three  under  any 
constellation.  Like  the  4x4  quasi-orthogonal  codes  presented  in  the  previous  section,  this  code  is  also  a  linear 
dispersion  code  constructed  from  HR  matrices.  This  8x4  unit-rate  code  is 

C(*i>*2)=Mo*i  +  ]A^x2,Aix1  +)A5X2,A2XX  +  ]A6X2,Aixx  +j^7A2],  (44) 

v - v - '  " - v - '  " - v - '  ' - v - ' 

£i  £2  £3  £4 


where  x,  and  x2  are  two  real-valued  8x1  symbol  vectors,  and  Aj,  i  =  0, . . . ,  7,  are  8x8  matrices  from  any 
single  HR  family  of  size  eight  satisfying  (14).  This  code  is  motivated  by  the  structure  of  a  half-rate  4x4 
orthogonal  code  for  four  transmitters.  In  fact,  it  is  easy  to  verify  that  the  code  (44)  is  an  orthogonal  code  if  the 
8x1  complex  vector  x{  +  jx2  is  replaced  by  a  4  x  1  real  vector  and  At,  i  =  0, . . . ,  3  are  replaced  by  4  x  4  HR 
matrices.  This  code  has  a  very  simple  and  appealing  structure. 

The  above  code,  also  referred  to  as  non-orthogonal  HR  code,  is  a  special  form  of  a  more  general  HR  code 
introduced  in  [5],  and  is  also  a  special  form  of  the  linear  dispersion  codes  introduced  in  [29,34]. 
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One  specific  form  (example)  of  the  non-orthogonal  HR  code  follows  from  (44)  with  the  HR  family 


given  in  (7): 

/*i(l) 

+  j*2(6) 

xi(3)+jx2(7) 

xi(4)+jx2(8) 

xi(2)+jx2(5)  \ 

xi  (2) 

-  Jy2(5) 

Yi(4)-jx2(8) 

— xi(3)+jx2(7) 

-xi(l)+jx2(6) 

xi(3) 

jY2(B) 

-xi(l)  +jx2(5) 

xi(2)+jx2(6) 

-Xi(4)— jx2(7) 

xi(4) 

+  )x2  (7) 

-xi(2)  —  jx2(6) 

-xi(l)+jx2(5) 

*i(3)  -jx2(8) 

xi  (5) 

+  j*2(2) 

X!(7)-jx2(3) 

xi(8)  -  jx2(4) 

— xi(6)  —  jx2(l) 

xi(6) 

-j*2(l) 

xi(8)+jx2(4) 

— xi(7)  —  jx2(3) 

xi(5)  —  jx2(2) 

xi  (7) 

-j*2(4) 

-xi(5)  -jx2(l) 

xi(6)  —  jx2(2) 

xi(8)  +  jx2(3) 

vxj(8) 

+  jx2(3) 

-xi(6)+jx2(2) 

-xi(5)  -  jx2(l) 

-xi(7)+jx2(4) } 

(45) 


where  xfk)  is  the  kt h  element  of  the  real  symbol  vector  xt. 

From  the  theorem  shown  next,  the  code  (44)  is  guaranteed  to  have  diversity  three  even  when  all  symbols  are 
independently  selected  from  any  constellation  (regardless  of  the  design  of  the  constellation).  To  our 
knowledge,  for  four-transmitters  systems,  the  code  (44)  is  the  only  known  unit-rate  linear  dispersion  code  that 
is  guaranteed  to  be  of  diversity  (at  least)  three  under  any  given  constellation.  This  is  a  useful  property  in 
practice  since  any  symbol  constellation  can  be  applied  to  this  code  while  the  diversity  is  guaranteed  to  be  no 
less  than  three. 

Theorem  4.1.  Given  any  xx  +  jx2  A 0,  the  code  C{xx,x2)  defined  in  (44)  has  a  rank  no  less  than  three  regardless  of 
the  constellation  from  which  all  symbols  are  independently  selected. 

The  rest  of  this  section  is  to  prove  Theorem  4.1.  Since  the  proof  is  quite  lengthy,  we  divide  the  proof  into 
several  sections  as  explained  next. 

4.1.  Outline  of  the  proof  of  Theorem  4.1 

The  proof  consists  of  a  sequence  of  lemmas,  and  these  lemmas  are  indexed  as  follows: 

(  Lemma  4. 1  4=  Lemma  4.1.1 
Lemma  4.2  4=  Lemma  4.1.1 

'  Lemma  4. 3. 1.1  4=  (47), 

(14). 


Theorem  4.1  < 


Lemma  4.3  4=  Lemma  4.3.1 


All  lemmas  are  stated  below.  The  proofs  are  given  in  the  subsequent  subsections  in  the  order  shown  above. 
Theorem  4.1  results  immediately  from  the  three  main  lemmas: 

Lemma  4.1  ( Proof  in  Section  4.3).  The  minimum  rank  of  (44)  is  no  larger  than  three. 

Lemma  4.2  {Proof  in  Section  4.4).  Among  any  three  column  vectors  in  (44),  at  least  two  of  them  are  independent. 

Lemma  4.3  {Proof  in  Section  4. 7).  Given  (14),  if  any  three  column  vectors  in  (44)  are  linearly  dependent,  they  are 
orthogonal  to  the  fourth  vector  in  (44). 

The  above  three  main  lemmas  are  based  on  the  following  supporting  lemmas: 

Lemma  4.1.1  {Proof  in  Section  4.2).  Given  distinct  i,  j,  m,  n  and  the  {non-zero)  real  vectors  Xi  and  x2, 
the  equation  AjX{  +  jAjx2  =  k{Amxx  +  j Anx2)  holds  if  and  only  if 

k=k\  +  )k2  =  ±j, 

(AiAmAj  +  An)x  2  =  0,  (46) 

X]  =  —k2Aj  Anx2, 

where  k\  and  £2  are  real  numbers.  There  is  always  a  non-zero  x2  that  satisfies  {AjA^Aj  +  An)x2  =  0. 
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Lemma  4.3.1  {Proof  in  Section  4.6).  Referring  to  (44),  v2  is  orthogonal  to  v3  and  v4  if  (a)  the  condition  (14)  holds 
and  (b)  there  exist  two  complex  scalars  k=k\  +  )kj  and  t=t\  +  j?2  with  \t\ 2  +  \k\2^0  and  two  real  vectors  xx  and 
x2  such  that 


Hi 


(47) 


Lemma  4.3. 1.1  {Proof  in  Section  4.5).  The  necessary  and  sufficient  condition  for  (47)  to  hold  is 


Mxx  =  0, 

x2  =  —{k\  +  t?fi  l{A^k2  +  A-jt2)T{Ao  —  A2k\  —  A^t^x^, 


(48) 


where 


M  =  A2k2  +  A3t2  +  {k2  +  t2)  1  {A4  —  A§k\  —  Ajti){A(,k2  +  A2t2)T{Ao  —  A2k\  —  A2tx).  (49) 

One  important  property  of  M  is 

MtM  =  cMI  +  N,  (50) 


where  cm> 0,  N  is  symmetric  and  orthogonal,  and  the  eigenvalues  of  N  are  ±cm  of  differing  signs.  The  exact 
content  of  cm  is  given  in  (68)  and  the  content  of  N  is  given  in  (68)  and  (65).  Furthermore,  M  is  singular  if  and 
only  if 


k2~\~t2  —  k\  +  t\  +  1, 
t\ti  +  k\k2  =  0. 


(51) 


4.2.  The  proof  of  Lemma  4.1.1 

Given  the  complex-valued  equation  A,-x x  +  jAjx2  =  k{Amxx  +  j Anx2),  there  are  two  corresponding  real¬ 


valued  equations  (real  and  imaginary  parts): 

{A,-  -  k\Am)xx  +  k2Anx2  =  0,  (52) 

{Aj  -  k\An)x2  -  k2Amxl  =  0.  (53) 

Because  At  —  k\Am  is  an  invertible  matrix  for  distinct  i  and  m,  has  to  be  non-zero.  Otherwise,  there  is  no 
non-zero  solution  for  xx  or  x2.  From  (53),  we  have 

=  k2xATm{Aj  -  k\An)x2.  (54) 

Taking  (54)  into  (52)  leads  to 

Zx2  =  0,  (55) 

where  Z  =  {At  —  k\Am)k2l Ajn{Aj  —  k\An)  +  k2An.  It  is  clear  that  Z must  be  singular,  or  otherwise  there  is  no 
non-zero  solution  for  x2. 

With  Properties  2  and  4  of  the  HR  matrices,  it  is  easy  to  verify  that 

ZTZ  =  [k\  +  (1  +  k2)2k22  +  2  k2]h  +  2  AjAmAjAn.  (56) 


It  is  known  that  Z  is  singular  if  and  only  if  Z  has  at  least  one  zero  eigenvalue.  Therefore,  based  on  (56)  and 
Property  5,  the  matrix  Z  is  singular  if  and  only  if  k\  +  (1  +  k2)2k22  +  2k\  =  ±2.  This  equation  is  equivalent  to 
any  of  the  following  equations: 

k\  +  (1  +  k2y  +  2  k\k\  =  ±2kj, 
kl  +  2  (k2  +  1  )A~2  +  (1  +  k2)2  =  i4^2. 

(1  +k]  +  k\)2  =  ±4k\. 
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From  the  above,  we  see  that  the  minus  sign  leads  to  no  real  solution.  After  dropping  the  minus  sign,  the  above 
is  equivalent  to  any  of  the  following: 

{\+k]+k\-  lk2)(\  +  k]  +  k\  +  2k2 )  =  0, 

(kj  +  (1  —  k2)2)(k2  +  (1  4-  k2)2)  =  0, 

k\  =  0  and  k2  =  ±  1 .  (57) 

Therefore,  k  =  ±j.  Taking  this  back  into  (54)  and  (55)  yields  the  sufficient  and  necessary  conditions  on  X2  and 
xx  as  shown  in  the  lemma. 

4.3.  The  proof  of  Lemma  4.1 

It  follows  from  Lemma  4.1.1  that  there  are  always  non-zero  xx  and  x2  such  that  any  two  columns  from  (44) 
are  linearly  dependent  of  each  other.  For  example,  if 

*1  =  [-l,  1,-1, 1,-1, -1,-1, -1]T, 

*2  =  [-l,  U, -1,-1, -1,1,  If 

then  the  codeword  given  in  (44)  has  a  rank  no  more  than  three. 


4.4.  The  proof  of  Lemma  4.2 


Based  on  Lemma  4.1.1  we  can  now  prove  that  any  three  columns  from  (44)  have  a  rank  larger  than  one. 
Suppose  that  Atxx  +  jAjx2  depends  on  each  of  Amxx  +  )Anx2  an<3  Atxx  +  }Arx2.  From  Lemma  4.1.1,  we  have 

(  x}  =  —k2AjAnx2, 

l  =  (58) 


which  implies  Aj(An  ±  Ar)x2  =  0.  However,  Aj(An  ±  Ar )  is  an  orthogonal  matrix,  which  means  that  x2  =  0 
and  hence  xx  =  0.  Therefore,  Aixx  +  ]Ajx2  cannot  depend  on  each  of  Amxx  +)Anx2  and  Atxx  +j^4rx2- 


4.5.  The  proof  of  Lemma  4. 3. 1.1 


Given  the  complex- valued  equation  (47),  we  equivalently  have  the  following  two  real- valued  equations 
(i.e.,  the  real  and  imaginary  parts): 

A0xx  =  A2xxk\  +  ^3*1/1  —  A()X2k2  —  A2x2t2, 

A4X2  =  A2xxk2  +  A2xxt2  +  A(,x2kx  +  A2x2tx 
or  equivalently 

(A0  -  A2kx  -  A3tx)xx  =  -(A6k2  +  A7t2)x2,  (59) 

(A4  —  Afrkx  —  Ajtx)x2  =  ( A2k2  +  A2t2)xx.  (60) 

Recalling  Property  4,  (Am  —  Ankx  -  Aitx )  is  always  non-singular  for  distinct  m,  n  and  /,  and  ( Amk2  +  Ant2)  is 
non-singular  for  distinct  m  and  n  unless  k\  +  t\  =  0.  So  from  (59)  and  (60),  k2  +  t2^0  unless  both  xx  and  x2 
are  equal  to  zero.  Also  from  (59)  and  (60),  xx  =  0  if  and  only  if  x2  =  0. 

From  (59),  we  have 

x2  =  ~{k\  +  t\)~l{A6k2  +  A7t2)T(A0  -  A2kx  -  A3tx)xx.  (61) 

Also,  from  (60),  we  have  an  equivalent  form  of  x2: 

x2  =  (1  +kx  +  t\)  l(A4  —  A(,kx  —  A7tx)T(A2k2  +  A3t2)xx.  (62) 
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Using  (61)  in  (60)  yields  Mxx  =  0  where 

M  =  M\  +  M2  (63) 

with  M i  =  A2k2  +  A3t2  and 

Mi  =  (k\  +  t2)  l(A4  —  A(,k\  —  Ajtx)(A^k2  +  Ajt2)T(Ao  —  A2k\  —  A3tx). 

Clearly,  xx  £  0  if  and  only  if  M  is  singular.  Summarizing  the  above,  we  have  that  (47)  holds  if  and  only  if  M  is 
singular,  xx  satisfies  Mxx  —  0,  and  x2  satisfies  (61). 

To  reveal  a  property  of  MTM,  we  now  apply  Property  4  to  obtain  that  MXM\  =  (k\  +  t2)I%  and 
M2M\  =  (k\  +  1  (1  +  k\  +  tl)2I%.  Furthermore, 

MjM2  =  (k\  +  t2)  l(Ajk2  +  Aj t2)(A4Ajk2  +  A4AJ t2  —  kxk2I%  —  A^A^kxt2  —  A2Ajtxk2  —  txt2I%) 
x (A0  -  A2kx  -A3tx) 

=  Mc  12,1  +  Mc  12,2,  (64) 

where 

Mc  12,1  =  ( k\  +  t\ )  l(Ajk2  +  Aj t2)  [A4AJ k2  +  A4A^t2  +  ^6^7  (— +  ^2^i)](^o  —  ^2^1  —  A3tx),  (65) 

' - - - " 


Mc  12,2  —  —  (k\  +  ^)  !(^42^2  +  ^3 t2)(txt2  +  kxk2)(Ao  —  A2kx  —  A3tx) 

=  —  (k2  +  t2)  \txt2  +  kxk2)(Ajk2Ao  —  kxk2Ig  —  AjA3txk2  +  Ajt2Ao  —  Ajt2A2kx  —  txt2Is).  (66) 
Using  Property  2  and  (66)  yields 

Mc  12,2  +  M^X22  =  2(k\  +  ^)  H^U2  +  k2ki)~I%.  (67) 

Therefore, 

MJM  =  M]M\  +  M2  +  Mc12,i  +  Mci2,2  +  Mj12jl  +  Mj12j2 

=  [(^2  +  4)  +  (^2  +  ^2)  *0  +  ^1  +  ^i)2  +  2{k\  +  Z^)  l(kxk2  +  ?i^)2]78  +  Mc12,i  +  MtcX2,x  .  (68) 


It  is  straightforward  to  verify  that 
TV  =  Mc  12,1  +  MtcX2X 

=  2(k\  +  ?2)_1  —  /r2^2^^?i^3  —  ^3  h^Mk\A2] .  (69) 

- - - v - 1 - - - ' 

Wo 

We  note  that  =  caIs  where  ca  =  t2  +  k\  +  (k2tx  —  t2kx)2.  It  is  also  straightforward  to  verify  that 

W0WT0=c2Ah.  (70) 

We  now  need  to  prove  Wq£  ±  c^/g.  Suppose  Wq  =  ±c^/g.  Then, 


AJ^AmAq  +  AJ^AmAo  —  k2AlAMt\A3  —  Alt2AMk\A2  =  ±caI% 


which  can  be  rewritten  as 

k2Aj(Ao  —  A3tx)  =  ±/(/1m)T  —  t2Aj(Ao  —  kxA2),  (71) 

where  f(AM )  =  AjA6k2  +  AjA7t2  +  A^A^-kAi  +  k2tx). 

Note  that  f(AM )  is  antisymmetric,  and  f(AM)Tf(AM)  =  cAh-  For  Ak  where  k  £(4,6,1),  AJAm  =  f(AM)Aj.. 
From  (71),  we  have 

{k\  +  t22)A0  =  ( A2k2  +  A2t2){±f{AM)  -  A]A2(k2tl  -  klt2)}.  (72) 
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We  now  multiply  (72)  by  Aj  to  yield 

— (k\  +  t2)AjAo  =  —  (Aj  A2k2  +  ^7){±/T(d^)  —  AjA2(k2tx  ~  kxt2)} 

=  ±  [k2Aj A2AJ Ae  +  k2t2AjA2AjA2  +  k2{k2t\  — 

±  +  (— 1)^2 (^2^1  —  k\t2)I%  +  AjA2(k2t\  —  kxt2)t2-  (73) 


Applying  Properties  2  and  5,  we  add  Eq.  (73)  to  its  transposed  version  to  yield 


\ 


k2  ±(A:2^4j^42^4^46  +  ^2^43^2^44^47  +  (kit\  —  k\t2)Aif)  A2AJAj)  +  (—l)(k2ti  —  kxt2)Is 


Fo 


=  0. 


/ 


(74) 


Recall  kl  +  tl^O-  From  the  definition  of  F0  shown  in  (74),  it  is  easy  to  verify  that  FqFJ  = 
[k\  +  t\  4-  (k2t\  ~  t\k2)2]h  :A(k2t\  —  t\ki)2-  Therefore,  (74)  implies  k2  =  0.  Similarly,  we  can  multiply  (72)  by 
Aj,  and  then  follow  the  same  analysis  as  shown  above  to  conclude  that  (2  =  0.  Therefore,  we  have  proven  by 
contradiction  that  Wo  A1  ±  c^/g. 

Since  MJM  =  cm! 8  +  N  where  N  has  the  eigenvalues  ±2^  +  A )  Xca,  M  is  singular  if  and  only  if 

cm  =  2(k\  T  ^2)  1  ca  (75) 

which  is  equivalent  to 

[(k2  +  tlf  ~  (*1  +  t\  +  l)2]2  +  4(txt2  +  kik2)2  =  0 

and  hence  (51). 


4.6.  Proof  of  Lemma  4.3.1 

Given  the  definitions  of  v2  and  v3  shown  in  (44),  it  follows  that 
v2v3  =  QcJaJA2x1  +  xJaJA6x2)  +  (xJaJA6x2  -  xjA]A2xl)j 

=  (xJaJA6x 2  -  xjAjA2xl)j.  (76) 

Using  (61)  and  Property  3,  the  first  term  (ignoring  j)  in  (76)  becomes 

xJaJA6x2  =  -  (k22  +  t2)~^ x^ Aj Ao(A^k2  +  ^7^)T(-4o  -  A2kx  -  ^i)*i 

=  —  (k2  +  t2)-lxJ(k2AjA0  —  k\k2A\A2  —  k2t\A\A3  +  t2AjAoAjAo 
—  t2kiAjA6AjA2  —  t2tiAjA6AjA3)x1 

=  ~(k\  +  t2)~lxj (t2AjA6AjA0  -  t2klAjA6AjA2  -  t2t\ATxA6AT1A3)xl.  (77) 

Taking  (62)  into  the  second  term  in  (76),  we  can  similarly  show  that 

x^2^5^2  =  (1  +  k\  +  ^)_1v^2^5(^4  -  A6kx  -  A1tx)T(A2k2  +  A3t2)xx 

=  (1  +  k\  +  t2)~]t2xT(AjA5AjA3  -  kxAjA5AjA3  -  txAjA5AjA3)xx.  (78) 

Based  on  (14)  and  Property  2,  we  have  AjA6AjAo  =  AjA2AjA3  = —AJA5AJA3.  Also,  recall 
1  +  k2  +  t\  =  k\  +  t2.  Therefore, 

—  (k-2  T  ^2 )  ^J(~LAjAoAjAo  +  t2kxAjA(,A2A2  +  t2 liAj AoAjA3 
—  ?2^42^5^4^3  +  kx  t2Aj  A3Aj  A3  +  tx  t2A^ A3Aj A3)x  J 
=  (k\  +  t2rlt2xTx(kxATxA6AT1A2  +  txAjA6AjA3  +  kxAT2A5AT6A3  +  txAjA5AjA3)xx  j.  (79) 

V - V - ' 

mn 
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In  order  to  prove  i§v3  =  0,  it  remains  to  prove  x^M^xx  =  0.  By  Lemma  4.3.1. 1,  Mx{  =  0  and 
MJ M  =  cmH  +  N.  From  the  property  of  N,  it  follows  that  range(c^/8  +  N)  and  range(cMfs  —  N)  are 
orthogonal  complement  of  each  other.  Therefore,  the  solution  space  of  Mxx  =  0  is  given  by  the  range  of 
Mq=cmI 8  —  N,  i.e.,  Xi  =  M0  v  for  any  real  vector  y.  Then,  xf  MNx}  =  0  if  and  only  if  v  =  0  for 

any  y.  The  proof  of  y7  M^MnMq  v  =  0  is  straightforward  but  very  lengthy,  the  details  of  which  are  given  in 
the  Appendix. 

To  prove  i§y4  =  0,  we  need  to  exchange  A2  with  A 3  and  A e  with  Aj  in  the  proof  for  yfy3  =  0.  With 
Property  6,  (14)  still  holds  after  the  double  exchanges.  So,  the  proof  of  =  0  is  basically  identical  to  the 
proof  of  t§y3  =  0. 

4.7.  Proof  of  Lemma  4.3 

Each  vector  in  [yl,y2,y3,y4]  depends  on  two  out  of  eight  HR  matrices  satisfying  (14).  Exchanging  any  two 
vectors  in  [v\,v2,v3,V4\  is  equivalent  to  exchanging  two  pairs  of  HR  matrices.  With  Property  6,  condition  (14) 
continues  to  hold  under  any  even  number  of  exchanges  of  HR  matrices.  Therefore,  following  the  same  proof 
as  for  Lemma  4.3.1,  if  any  three  vectors  in  \y\,y2,y3,y4\  are  linearly  dependent  of  each  other,  they  must  be 
orthogonal  to  the  fourth  vector. 


5.  Further  remarks  on  the  non-orthogonal  HR  codes 

5.1.  Full  diversity  non-orthogonal  HR  codes 

Like  the  quasi-orthogonal  codes,  the  non-orthogonal  HR  codes  can  also  be  made  full  diversity  by 
introducing  proper  diversity  in  symbol  constellations.  While  the  best  method  to  achieve  full  diversity  of  the 
non-orthogonal  HR  codes  is  still  an  open  problem,  we  give  one  method  here  to  achieve  full  diversity.  Consider 
the  codeword  as  in  (45).  Let  each  complex  element  of  x,^)  +  j x2(i),  i=  1,2, 3,4,  be  from  the  constellation  set 
of  exp[j7i(m/4  +  1/6)],  m  =  1, 3, 5, 7,  and  each  complex  element  of  x,(i)  +  jx2(z),  i  =  5, 6, 7, 8,  be  from  another 
constellation  set  exp(j7im/4),  m  =  1, 3, 5, 7.  Through  exhaustive  search,  it  has  been  verified  that  the  rank  of  the 
four  column  code  matrix  is  always  four. 


5.2.  Simulation 

To  illustrate  the  performance  of  the  non-orthogonal  HR  codes,  we  show  a  simulation  example.  For  a 
system  of  single  receiver  and  four  transmitters,  each  block  of  received  data  can  be  expressed  as 

y  =  C(x)  h  +  n ,  (80) 

where  we  assume 

•  The  entries  of  the  fading  vector  h  are  i.i.d.  Gaussian  distributed  with  unit  variance. 

•  The  entries  of  the  noise  vector  n  are  i.i.d.  Gaussian  distributed  of  variance  x  4  x  io_01SNR  where  SNR  is 
the  dB  value  of  the  ratio  of  the  transmitted  power  over  the  noise  variance.  Here,  is  dependent  on 
modulations.  For  example,  u2m  =  2  for  4-QAM  (four  symbols  on  the  corners  of  a  square  of  side  equal  to  2), 
o2m  =  1  for  QPSK  (four  symbols  uniformly  spaced  on  a  circle  of  unit  radius),  and  o2m  =  10  for  16-QAM. 
The  factor  4  is  due  to  four  transmitters. 

•  The  code  matrix  C(x)  is  given  by  (45). 

An  alternative  form  of  (80)  is  as  follows: 

(  9*00  \  /  X\  \  (  9*(«)  \ 

(81) 
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Fig.  2.  Averaged  bit  error  rate  versus  SNR  for  four  space-time  block  codes.  The  data  rate  for  all  codes  is  2bits/s/Hz.  A  system  of  four 
transmitters  and  one  receiver  is  considered. 


where 


//  = 


E;=o  i) 


ELo^'+4  aw\ 
Ew^SKW/’ 


(82) 


We  used  (81)  with  a  sphere  decoding  algorithm  [35-37]  to  detect  v.  For  each  realization  of  x,  we  chose  an 
independent  realization  of  h  and  n. 

In  our  simulation,  we  compared  the  four  different  codes:  the  quasi-orthogonal  code  [11],  and  the  full  rank 
quasi- orthogonal  code  with  the  constellation  rotation  given  in  [14,16],  the  non-orthogonal  code  (44),  and 
the  half-rate  complex  orthogonal  code  [9].  To  ensure  the  same  bit  rate,  we  used  QPSK  for  the  first  two  codes, 
4-QAM  for  the  third  code,  and  16-QAM  for  the  fourth  code. 

The  performances  of  the  four  different  codes  are  compared  in  Fig.  2.  The  non-orthogonal  HR  code 
(referred  to  as  diversity  3  code  in  the  figure)  shows  a  good  performance  in  a  medium  range  of  SNR,  i.e.,  better 
than  the  original  quasi-orthogonal  code  and  even  the  half-rate  orthogonal  code.  The  half-rate  orthogonal  code 
performs  well  at  very  high  SNR  because  of  its  full  diversity.  But  the  full  diversity  quasi-orthogonal  code 
performs  the  best  among  the  four  codes  compared. 

It  is  important  to  remember  that  the  code  (44)  has  diversity  no  less  than  three  for  any  constellation  while  the 
full  diversity  quasi-orthogonal  code  needs  to  be  readjusted  for  each  different  constellation. 


6.  Conclusion 

In  this  paper,  we  have  investigated  STBC  that  have  strong  connections  with  the  HR  families  of  matrices. 
The  key  contributions  are  Theorems  3.1  and  4.1.  Theorem  3.1  states  that  the  Type  I  family  of  all  published  as 
well  as  unpublished  4x4  unit-rate  quasi-orthogonal  codes  are  simply  variations  from  two  independent  codes 
shown  in  (16)  via  (15).  Theorem  4.1  states  that  the  unit-rate  code  (44)  for  four  transmitters  has  a  rank  no  less 
than  three  under  any  given  constellation.  To  our  knowledge,  the  code  (44)  is  the  only  known  unit-rate  linear 
dispersion  code  that  is  guaranteed  to  have  diversity  (at  least)  three  for  four  transmitters.  This  is  a  useful 
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advantage  since  it  could  reduce  the  physical  layer  complexity  associated  with  constellation  constraint.  It 
remains  a  challenge  to  discover  whether  or  not  there  exists  a  linear  dispersion  code  that  guarantees  a  higher 
diversity  than  the  code  (44)  for  four  transmitters  over  all  possible  constellations.  It  is  our  hope  that  the 
in-depth  analysis  shown  in  this  paper  will  motivate  and  help  this  pursuit. 
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Appendix  A.  Proof  of  v  =  0 

The  proof  is  relatively  lengthy.  We  will  repeatedly  apply  Properties  1-7  of  the  HR  matrices  as  well  as 
condition  (14).  To  help  the  presentation  of  the  proof,  we  will  use  the  sum  table  to  be  introduced  next. 


A.l.  Introduction  of  the  sum  table 

For  example,  to  describe  the  sum  ai AJA2AJA4  +  0C2AJAgAJA3,  we  use  the  following  sum  table: 


coefficients 

products  of  matrices 

Ot\ 

[1234] 

OL2 

[5673] 

(83) 


where  each  group  of  four  integers  from  (0, 1,2, 3,4, 5,6, 7)  corresponds  to  a  product  of  four 
matrices,  i.e., 

±[ijmn]=  ±  AjAjAjnAn.  (84) 

To  describe  filAjA4(aiAjA2AjA4  +  a2AjA^Aj A3)  as  another  example,  we  will  use 


coefficients 

Products  of  matrices 

0i&i 

/?lCt2 

[24]  [1234] 

[24]  [5673] 

From  Properties  1(a)  and  2(a),  it  is  easy  to  verify  that  [241234]  =  [13].  Because  of  AjA^AjAo  = 
AJA2AJA3  =  AjA4AjA3,  we  can  write  [245673]  =  [245367]  =  [1670][67]  =  [01].  Therefore,  we  can  also 
express 

^1AjA4(a1AjA2AjA4  +  oc2AJA6AJA3) 


as 


coefficients 

Products  of  matrices  and  simplification 

PlOL  1 

PlO/,2 

[24]  [1234]  =  [13] 

[24]  [5673]  =  [245367]  =  [1670]  [67]  =  [01] 
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A.2.  Main  body  of  the  proof 

From  the  definition  of  Mo,  it  is  straightforward  to  verify  that 
M0  =  cmI  8  —  N 

=  2{k22  +  t2)I%  ~  2k2(k2  +  1 2)  lAj[A4Aj  k2  +  ^4^7  ^2  +  A^Aj(— k\t2  +  k2t\j\(Ao  —  A^tx) 
—  2{k2  +  t2)  +  A4A2t2  +  AfiAj(— k\t2  +  k2tx)\(Ao  —  k\A2). 

With  c\=k\  +  t\  and  c2=kxt2  +  k2tx,  we  can  write 


(87) 


M0  = 


coefficients 

products  of  matrices 

2Cl 

h 

— 2k2C±l 

[2460] 

-2k2t2c^1 

[2470] 

[2670] 

2k\t\c{1 

[2463] 

2k2t2t\c{x 

[2473] 

2k2tic2c f1 

[2673] 

-2k2t2c{x 

[3460] 

-2t2t2c{x 

[3470] 

—  2c2^2Ci  1 

[3670] 

2k\k2t2c[l 

[3462]  =  -[2463] 

2kit2t2c^1 

[3472]  =  -[2473] 

2kic2t2c^1 

[3672]  =  -[2673] 

(88) 


Therefore, 


M0  = 


coefficients 

products  of  matrices 

2C! 

h 

-2  k^1 

[2460] 

-2k2t2c{x 

[2470] 

—2  k2c2c{x 

[2670] 

2k2tiCi 1  —  skik^c^1  =  2c2k2Ci1 

[2463] 

2k2t2t\c i1  —  2k\t2t2c[l  =  2c2t2c[1 

[2473] 

2k2t\c2c{ 1  -  2kic2t2c^1  =  2c2c^1 

[2673] 

—2k2t2cf1 

[3460] 

-2  t2t2c^ 

[3470] 

-2c2t2C]“1 

[3670] 

(89) 
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From  (79)  we  have 


M^y  —  k\AjA^AjA2  +  1\AJA^AJAj  +  /ciAjA^AjA^  +  tiAjA^AjAj 
=  k\[\612\  +  /i[1673]  +  ^[2563]  +  ti[2573]. 


(90) 


To  compute  the  product  MqMn  (where  M0  is  symmetric),  we  have 


M0  x  A;i  [1672]  = 


coefficients 

products  of  matrices  and  simplification 

group  index  in  (95) 

2k\C\ 

h  [1672]  =  [1672] 

1 

—ZkikzC^1 

[2460]  [1672]  =  [4017] 

2 

-2kik2t2c[1 

[2470]  [1672]  =  -[4016] 

3 

—2kik2C2c[1 

[2670]  [1672]  =  -[01] 

4 

2kic2k2c^1 

[2463]  [1672]  =  [4317]  =  [6025] 

because  [1670]  =  [5243]  [5620]  =  [1743] 

5 

2k\c2t2c^1 

[2473]  [1672]  =  -[4316]  =  [7025] 

because  [1670]  -  [5243]  [5270]  -  [1643] 

6 

2kiC2C^1 

[2673]  [1672]  =  -[31] 

7 

-2kik2t2c[l 

[3460]  [1672]  =  -[3460]  [5043]  =  [65] 

8 

-2k\t2t2c{x 

[3470]  [1672]  =  -[3470]  [5043]  -  -[57] 

9 

-2A:ic2t2c]"1 

[3670]  [1672]  =  -[3012] 

10 

(91) 


M0  x  ti  [1673]  - 


coefficients 

products  of  matrices  and  simplification 

group  index  in  (95) 

2iici 

Is  [1673]  =  [1673] 

11 

~2t\ klc^1 

[2460]  [1673]  =  -[2460]  [5240]  =  [56] 

8 

— 2t\k2t2c^^ 

[2470]  [1673]  =  -[2470]  [5240]  -  [75] 

9 

—2tik2c2c{ 1 

[2670]  [1673]  w  -[2013]  =  [3012] 

10 

2tic2k2c^1 

[2463]  [1673]  =  [2417]  =  -[6053] 

because  [1670]  =  [5243]  [5630]  =  [1247] 

12 

2t\c2t2  cjf1 

[2473]  [1673]  =  -[2416]  =  -[7053] 

because  [1670]  =  [5243]  <K>  [5703]  =  [1624] 

13 

2t\C2C^1 

[2673]  [1673]  =  [21] 

14 

-2tik2t2c{1 

[3460]  [1673]  =  [4017] 

2 

-2t\t2t2  c]-1 

[3470]  [1673]  =  -[4016] 

3 

2t\c2t2c^^ 

[3670]  [1673]  =  -[01] 

4 

(92) 
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coefficients 

products  of  matrices  and  simplification 

group  index  in  (95) 

2k\C\ 

7s  [2563]  =  [2563] 

2 

— 2&i&2cr1 

[2460]  [2563]  =  -[4053] 

1 

-2/ciA^c^1 

[2470]  [2563]  =  [2470]  [1407]  =  -[12] 

14 

— 2/ci/c2c2cf 1 

[2670]  [2563]  -  [7053] 

13 

2kic2k2c±1 

[2463]  [2563]  =  [45] 

15 

2k\c2t2c^1 

[2473]  [2563]  =  -[4756]  =  -[3012] 

because  [1670]  =  [5243]  <^>  [1230]  =  [5647] 

10 

2kiC2c{1 

[2673]  [2563]  =  -[75] 

9 

-2k\k2t2c[l 

[3460]  [2563]  =  -[4025] 

11 

—2kit2t2c^ 1 

[3470]  [2563]  =  [3470]  [1407]  =  [13] 

7 

-2/ciC2t2c]_1 

[3670]  [2563]  =  [7025] 

6 

(93) 


Mo  x  [2573]  = 


coefficients 

products  of  matrices  and  simplification 

group  index  in  (95) 

2t\C\ 

is  [2573]  =  [2573] 

3 

— 2tik%Ci1 

[2460]  [2573]  =  -[2460]  [1406]  =  -[21] 

14 

-2t\k2t2c{x 

[2470]  [2573]  =  -[4053] 

1 

-2ti/c2c2c]1 

[2670]  [2573]  =  -[6053] 

12 

2t\c2k2c [f1 

[2463]  [2573]  =  -[4657]  =  [3012] 

10 

2t\c2t2c{v 

[2473]  [2573]  =  [45] 

15 

2fic|c]f1 

[2673]  [2573]  =  [65] 

8 

— 2£iA;2£2c1  1 

[3460]  [2573]  =  -  [3460]  [1406]  =  [13] 

7 

—2  tit2t2c±l 

[3470]  [2573]  =  -[4025] 

11 

—2t\c2t2c^1 

[3670]  [2573]  =  -[6025] 

5 

(94) 


The  following  table  shows  how  the  common  terms  in  the  previous  four  tables  are  combined. 
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Group  index 

Final  Sum 

1 

2/ciCi[1672] 

2 

2/ciCi[2563] 

3 

2t±ci  [2573] 

4 

0 

5 

0 

6 

0 

7 

-2fcic|q1[31]  +  (-2fcit2t2cr1)(-[13])  +  (-2t1fc2t2q1)[13] 

8 

— 2kik2t2c^1[Qh]  +  —  2£i&2C^1[56]  +  2tiC2C^1[65] 

9 

2k]t2t‘2Cx  1 1 5  7 1  +  2(]  A:2/2e,  1 1 7  5 1  +  —2k\cr2c.x  1 1 7  5 1 

10 

0 

11 

2tiCi[1673] 

12 

0 

13 

0 

14 

2Gc|c^1[21]  +  2/cifc2f2cr1[12]  +  2t1k%ci1[21] 

15 

0 

Therefore, 


coefficients 

products  of  matrices 

2k\C\ 

[1672] 

2k\C\ 

[2563] 

2t\Ci 

[2573] 

2t\C\ 

[1673] 

2&ic2cjf1  +  2kit2Ci 1  —  2tik2t2c^1 

[13] 

2txc\c{1  +  2G&2V  —  2kik2t2c^1 

[65] 

2k\  c2c^  +  2/cit2c i1  —  2t\k2t2Ev1 

[57] 

— 2Gc2cj_1  —  2 tik^c^1  4-  2k\k2t2E[ 1 

[12] 

(96) 


We  now  need  to  multiply  each  term  of  MqMn,  i.e.,  (96),  by  M o  from  right.  Note  that 
vT(AjAiAjAtAjAr)v  =  0  and  vt(AJA;)v  =  0  for  any  vector  v  and  distinct  indices.  So,  we  will  use  [kijtqr]~ 0 
and  [ki\~ 0  to  denote  that  the  corresponding  terms  become  zero  after  being  multiplied  by  v  from  left  and  right. 
We  will  use 


E\=2k\c\cx  1  +  2k\t\cx  1  —  2t\k2t2Cx  \  (97) 

E2=2t\c\cxx  +  2t\k\c^1  —  2kxk2t2cx\  (98) 
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2fcici[1672]  x  M0  = 


coefficients 

products  of  matrices  and  simplification 

group  index 

2&ici2ci 

[1672]I8  =  [1672] 

1 

— 2k\C\2k2C^1 

[1672]  [2460]  =  [1740] 

2 

—2kiCi2k2t2Ci1 

[1672]  [2470]  -  -[1640] 

3 

—2kiCi2k2C2C^1 

[1672][2670f«  -[10]  ~  0 

2kici2c2k2c^1 

[1672]  [2463]  =  [1743] 

4 

2k\C\2c2t2c[l 

[1672]  [2473]  =  -[1643] 

5 

2kiC\2c2C ^ 

[1672][2673J&  -[13]  ~  0 

^kic^k^c^1 

[1672]  [3460]  m  [172340]  ~  0 

—2k\C\2t2t2c^  i 

[1672]  [3470]  =  -[162340]  ~  0 

— 2kiCi2c2t2Ci1 

[1672]  [3670]  =  -[1230] 

6 

(99) 


2fciCi[2563]  x  M0  = 


coefficients 

products  of  matrices  and  simplification 

group  index 

2k\C\2c\ 

[2563]  I8  =  [2563] 

2 

—2k\Ci2k2CiX 

[2563]  [2460]  =  -[5340] 

1 

— 2k\Ci2k2t2c[1 

[2563]  [2470]  =  -[563470]  ~  0 

— 2kici2k2c2c^1 

[2563]  [2670]  =  [5370] 

7 

2kici2c2k2c i1 

[2563]  [2463]  =  [54]  ~  0 

2kici2c2t2c^1 

[2563]  [2473]  =  -[5647] 

6 

2kiC\2c2Ci 1 

[2563]  [2673]  =  [57]  ~  0 

- 2kici2k2t2c ]1 

[2563]  [3460]  =  -[2540] 

8 

—  2^1012^2^2^  1 

[2563]  [3470]  -  [256470]  ~  0 

— 2k\  ci  2c2f  2  Ci1 

[2563]  [3670]  =  [2570] 

5 

(100) 
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2£ici[2573]  x  M0 


2tici[1673]  x  M0 


coefficients 

products  of  matrices  and  simplification 

group  index 

2£ici2ci 

[2573]J8  =  [2573] 

3 

— 2fiCi2/c2C^1 

[2573]  [2460]  =  -[573460]  -  0 

— 2iiCi2^2^2c]1 

[2573]  [2470]  -  -[5340] 

1 

-2tiCi2A:2C2c]1 

[2573]  [2670]  =  -[5360] 

9 

2fici2c2/c2c)1 

[2573]  [2463]  =  -[5746] 

6 

2£iCi2c2£2C^1 

[2573]  [2473]  =  [54]  ~  0 

2fici2c2cj 1 

[2573]  [2673]  =  [56]  ~  0 

- 2t1ci2k2t2Ci 1 

[2573]  [3460]  -  [257460]  -  0 

- 2tiCi2t2t2C i1 

[2573]  [3470]  =  -[2540] 

8 

—2t\C\2c2t2C\1 

[2573]  [3670]  =  -[2560] 

4 

coefficients 

products  of  matrices  simplification 

group  index 

2£ici2ci 

[1673]  Js  =  [1673] 

8 

— 2fiCi2fc|c^1 

[1673]  [2460]  =  [173240]  ~  0 

-2tiCi2k2t2C^1 

[1673]  [2470]  =  -[163240]  ~  0 

- 2t\C\2k2c2c i1 

[1673]  [2670]  =  -[1320] 

6 

2tic\2c2k2c Jf1 

[1673]  [2463]  -  [1724] 

9 

2fiCi  202^2^1  1 

[1673]  [2473]  =  -[1624] 

7 

2fici2c2c]1 

[1673]  [2673]  -  [12]  ~  0 

- 2tiCi2k2t2C j-1 

[1673]  [3460]  =  [1740] 

2 

—  2tiCi2f2^2ci  1 

[1673]  [3470]  =  -[1640] 

3 

-2^ici2c2t2c]“1 

[1673]  [3670]  =  -[10]  ~  0 

(101) 


(102) 
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coefficients 

products  of  matrices  and  simplification 

group  index 

Ei2ci 

[13 ]I8  =  [13]  ~  0 

-E12klcY 

[13]  [2460]  =  [132460]  ~  0 

-Ei2k2t2C^1 

[13]  [2470]  =  [132470]  ~  0 

-Ei2k2C2C^1 

[13]  [2670]  =  [132670]  ~  0 

E^c^Ci1 

[13]  [2463]  =  -[1246] 

7 

Ei2c2t2Ci 1 

[13]  [2473]  =  -[1247] 

9 

E^c^c^1 

[13]  [2673]  =  -[1267] 

1 

— £'l2/C2^2Ci  1 

[13]  [3460]  =  [1460] 

3 

—E^tfac^1 

[13]  [3470]  =  [1470] 

2 

—Ei2c2t2Ci 1 

[13]  [3670]  -  [1670] 

10 

coefficients 

products  of  matrices  and  simplification 

group  index 

Ei2ci 

[57 ]/„  =  [57]  ~  0 

-Ex2klc{x 

[57]  [2460]  =  [572460]  ~  0 

-Ei2k2t2E[l 

[57]  [2470]  =  [5240] 

8 

f- Ei2k2c2c^1 

[57]  [2670]  -  [5260] 

4 

E^c-z^c-l  1 

[57]  [2463]  =  [572463]  ~  0 

E12c2t2E[1 

[57]  [2473]  =  [5243] 

10 

E^c^c^1 

[57]  [2673]  -  [5263] 

2 

-Ei2k2t2E[l 

[57]  [3460]  =  [573460]  -  0 

—E^t^Ci 1 

[57]  [3470]  =  [5340] 

1 

— Ei2c2t2Ci  1 

[57]  [3670]  =  [5360] 

9 
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coefficients 

products  of  matrices  and  simplification 

group  index 

~E2  2ci 

[12]  Jg  =  [12]  ~  0 

E22kicY 

[12]  [2460]  =  [1460] 

3 

E22k2t2C1 1 

[12]  [2470]  =  [1470] 

2 

E^k^c^1 

[12]  [2670]  =  [1670] 

10 

—E'iZq^c-y  1 

[12]  [2463]  =  [1463] 

5 

—E22c2t2c-^ 1 

[12]  [2473]  =  [1473] 

4 

-Erfclci1 

[12]  [2673]  =  [1673] 

8 

E^k^c^1 

[12]  [3460]  =  [123460]  ~  0 

E22t2t2cx  1 

[12]  [3470]  =  [123470]  ~  0 

E22c2t2c^ 1 

[12]  [3670]  =  [123670]  ~  0 

(105) 


coefficients 

products  of  matrices 

group  index 

E22c\ 

[65]  J8  =  [65]  ~  0 

-E22klcY 

[65]  [2460]  =  -[5240] 

8 

-E22k2t2cx  1 

[65]  [2470]  =  [652470]  ~  0 

-E22k2c2c[1 

[65]  [2670]  =  [5270] 

5 

E22c2k2c^ 1 

[65]  [2463]  =  -[5243] 

10 

E22c2t2E[l 

[65]  [2473]  =  [652473]  ~  0 

E^clc^1 

[65]  [2673]  =  [5273] 

3 

—E^krf-zCi 1 

[65]  [3460]  =  -[5340] 

1 

E22t2t2c^ 1 

[65]  [3470]  =  [653470]  ~  0 

—£^2202^2^  1 

[65]  [3670]  =  [5370] 

7 

There  are  10  distinct  groups  of  common  terms  in  the  previous  eight  tables.  All  the  common  terms  cancel 
each  other  as  shown  by  the  following  tables. 
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The  1st  group: 


coefficients 

product  of  matrices 

2fcici2ci 

[1672] 

— 2fcici2&2dj"1 

-[5340]  =  -[1672] 

because  [1670]  =  [5243]  [1672]  =  -[5043] 

—2t\Ci2k2t2C[1 

-[5340] 

Ei2clc±l 

-[1267]  =  -[1672] 

—Ei2t2t2C1  1 

[5340] 

—E22k2t2C^  1 

-[5340] 

The  2nd  group: 


coefficients 

product  of  matrices 

— 2&iCi2&2c[f 1 

[1740] 

2k\C\2c\ 

[2563]  =  -[1740] 

because  [1670]  [5243]  [1470]  =  -[5263] 

-2tic12k2t2C^1 

[1740] 

— £’l2f2^2ci  1 

[1470]  =  -[1740] 

E^c^c^1 

[5263]  =  [1740] 

E22k2t2C^1 

[1470] 

The  3rd  group: 


coefficients 

product  of  matrices 

—2kiCi2k2t2C^1 

-[1640] 

2tici2ci 

[2573]  =  [1640] 

because  [1670]  =  [5243]  [1640]  =  -[5273] 

-2tiCi2t2t2C^1 

-[1640] 

- Ei2k2t2c i1 

[1460]  =  -[1640] 

E22klcY 

[1460] 

E22c%E[1 

[5273]  -  -[1640] 

(107) 


(108) 


(109) 
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The  4th  group  of  terms  in  the  sum  is 


coefficients 

product  of  matrices 

2kxcx2c2k2E[1 

[1743] 

— 2tici2c2t2c^1 

-[2560]  =  -[1743] 

because  [1670]  =  [5243]  [1473]  =  [5260] 

—E^k^Ci1 

[5260]  =  -[1743] 

—E^c^c-l  1 

[1473]  =  -[1743] 

The  5th  group: 


coefficients 

product  of  matrices 

2k1c12c2t2c^1 

-[1643] 

—2k1c12c2t2c^1 

[2570]  =  -[1643] 

because  [1670]  =  [5243]  [1643]  =  [5270] 

-Erfc^c^1 

[1463]  =  -[1643] 

—E^k-zc-z^  1 

[5270]  =  [1643] 

The  6th  group: 


coefficients 

product  of  matrices 

— 2kiCi2c2t2Ci 1 

-[1230] 

2/CiCi  202^2^  1 

-[5647]  =  -[1230] 

because  [1670]  =  [5243]  ^  [5674]  =  [1203] 

2t1c12c2k2c^1 

-[5746]  =  [1230] 

— 2tic±2k2c2Ci1 

-[1320]  =  [1230] 

The  7th  group: 


coefficients 

product  of  matrices 

-2/ciCi2fc2c2cf1 

[5370] 

2tici2c2t2<h  1 

-[1624]  =  -[5370] 

because  [1670]  =  [5243]  [5370]  =  [1246] 

£l2c2&2C]f1 

-[1246]  =  -[5370] 

—E22c2t2cx  1 

[5370] 
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(110) 


(111) 


(112) 


(113) 
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The  8th  group: 


coefficients 

product  of  matrices 

- 2k\C\2k2t2c if1 

-[2540] 

—2tiCi2t2t2C^1 

-[2540] 

2t1c12c1 

[1673]  =  [2540] 

because  [1670]  *  [5243]  44  [1673]  =  -[5240] 

-Ei2k2t2c^ 1 

[5240]  =  -[2540] 

-Ertclc^1 

[1673] 

-E22k22cY 

-[5240] 

The  9th  group: 


coefficients 

product  of  matrices 

-2fici2fc2c2cf1 

-[5360] 

2t1ci2c2k2c^1 

[1724]  =  -[5360] 

because  [1670]  =  [5243]  44  [1274]  =  [5603] 

E^crf^Y 1 

-[1247]  =  [5360] 

— £^i2c2f2c1 1 

[5360] 

The  10th  group: 


coefficients 

product  of  matrices 

— £^i2c2f2c1 1 
E\2c2t2c^ 1 
T^22fc2c2c]_1 
£l22c2fc2c1 1 

[1670] 

[5243]  =  [1670] 

[1670] 

-[5243]  =  -[1670] 

(114) 


(115) 


(116) 
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